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Penetration  of  Microwaves  Through  Dispersive  Concrete 
Using  a  Three-Dimensional  Finite-Difference  Time-Domain  Code 

S.  J.  Yakura  and  D.  Dietz 

Air  Force  Research  Laboratory,  Directed  Energy  Directorate 
Kirtland  AFB,  NM  87117-5776 

Abstract 

In  this  paper,  we  present  a  new  formulation  of  a  three-dimensional  Finite-Difference  Time-Domain  (FDTD) 
algorithm  that  is  applicable  for  studying  the  penetration  of  microwave  pulses  in  conductive  and  dispersive 
materials,  such  as  concrete.  We  investigate  the  ^ects  of  microwaves  penetrating  into  conductive  and  dispersive 
concrete  using  the  newly  formulated  algorithm  incorporated  into  a  three-dimensional  FDTD  code.  Because  of  the 
dispersive  properties  of  concrete,  both  conductivity  and  polarization  terms  are  taken  to  be  frequency-dependent.  We 
model  the  electric  polarization  terms  and  the  electric  conductivity  to  exhibit  linear  dispersive  behavior  via  first- 
order,  time-dependent,  electric  susceptibility  and  conductivity  functions.  Sample  one-dimensional  FDTD 
calculations  show  that  the  penetration  ofwide-band  and  ultrawide-band  microwave  pulses  must  be  simulated  using 
the  fully  dispersive  FDTD  algorithm  to  obtain  the  proper  response  from  dispersive  concrete.  On  the  other  hand,  the 
use  of  the  fully  dispersive  FDTD  algorithm  is  not  essential  in  simulating  the  penetration  of  narrow-band  microwave 
pulses  through  concrete  because  of  the  non-dispersive  nature  of  the  narrow-band  (single  frequency)  effect. 

1.  INTRODUCTION 

Determining  the  penetration  through  concrete  of  microwaves  using  numerical  modeling  has  received  much 
attention  in  recent  years.  Measured  ground  penetration  radar  data  in  the  rmge  of  600  MHz  to  3  GHz  has  shown 
concrete  to  exhibit  dispersive  characteristics  in  both  electrical  permittivity  and  conductivity.  Based  on  the 
experimental  data  now  available  from  many  experiments  that  have  been  carried  out  in  last  few  years,  we  have 
sufficient  information  to  simulate  the  detailed  behavior  of  microwave  pulses  that  penetrate  through  dispersive 

concrete  materials  [1-3].  .  .  „ 

With  the  advent  of  present  day  computers  which  provide  very  fast  execution  times  and  great  quantities  of 
computer  memory,  we  are  at  the  point  where  we  have  enough  computational  power  to  solve  Maxwell’s  equations 
directly  for  a  large-scale  problem.  Among  recently  investigated  numerical  techniques  that  show  great  promise  in 
achieving  this  goal  is  the  well-known  finite-difference  time-domain  (FDTD)  method  [4].  It  is  based  on  using  a 
simple  differencing  scheme  in  both  time  and  space  to  calculate  the  transient  behavior  of  electromagnetic  field 
quantities.  Because  of  the  simplicity  of  the  FDTD  method,  recent  researchers  have  focused  their  a^ntion  on  the 
numerical  evaluation  of  the  linear  convolution  integral  terms  which  appear  in  one  of  Maxwell’s  equations  (Ampere  s 
Law).  By  properly  evaluating  these  terms,  many  people  have  successfully  modeled  the  response  of  linear  dispersive 

effects  [5-10].  ^  • 

In  this  paper  we  consider  the  dispersive  materials  that  exhibit  the  linear  polarization  through  the  first-order,  time- 
dependent,  electric  susceptibility  function.  For  such  materials  the  relationship  between  D[t;x)  and  E(t;x)  can  be 
expressed  as  [1 1] 

^t;x)  =  €„€.  E(t;x)+e„  ^°\E(T;x)X^^Ht-T)dT ,  0-1) 

P-1  ^ 

where  €«  is  the  electric  permittivity  of  free  space,  e«  is  the  medium  permittivity  at  infinite  frequency,  and  (t-r) 
is  the  /Xh  term  of  the  collection  consisting  of  Pmax  first-order,  time-dependent,  elecnic  susceptibility  functions,  where 
Pmax  is  the  maximum  number  of  terms  which  we  choose  to  consider  for  a  particular  formulation  of  Eq.  (1.1).  ^ 

Also,  we  consider  that  the  current  density  vector,  J(t;x),  which  appears  in  one  of  Maxwell’s  equations  (Ampere’s 
Law),  to  be  dispersive.  To  include  the  dispersive  characteristics  of  the  current  density  vector,  we  first  divided  the 
current  density  vector  into  non-dispersive  and  dispersive  terms.  The  non-dispersive  term  is  expressed  in  the  usual 
form  as  the  product  of  the  constant  conductivity  coefficient  and  the  electric  field  vector,  E(t,x).  The  dispersive  term 
is  expressed  in  the  linear  convolution  integral  of  E(r;x)  and  the  first-order  time-dependent  conductivity  function. 
Hence,  the  expression  that  represents  the  current  density  vector,  J(t;x),  is  given  in  the  following  form: 


K  max  iCmax 

l(t;x)  =  E(t;x)'^CT*J^>  +  ^  jE(T;x)a<J^(t-T)dT  , 
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where  and  <yJ^^(t-T)  are  the  Kth  terms  of  the  collection  consisting,  respectively,  of  Kn^vc  constant  conductivity 
values  and  Kf„ax  first-order,  time-dependent,  conductivity  functions,  where  is  the  maximum  number  of  terms 
which  we  choose  to  consider  for  a  particular  formulation  of  Eq.  (1.2). 

By  including  the  dispersive  contributions  in  both  polarization  and  current  density  vectors  as  shown  in  Eqs.  (1.1) 
and  (1.2) ,  we  can  evaluate  the  linear  convolution  integrals  that  appear  in  Maxwell’s  equations  (Ampere’s  Law)  and 
come  up  with  a  second  order  accurate  FDTD  algorithm  to  simulate  the  characteristic  behavior  of  a  wave  propagating 
into  a  linear  dispersive  material.  We  investigate,  in  particular,  the  case  in  which  the  first-order,  time-dependent, 
electric  susceptibility  and  conductivity  functions  are  expressed  in  the  following  complex  exponential  forms  that 
contain  complex  constant  coefficients: 

X<J>(t)  =  Re{a^,exp(-'/;,t)}U(t),  (1.3) 

and 

=  Re  f  pj^exp  (  -6tt )}  U(t) .  (1.4) 

where  Ref  }  is  used  to  represent  the  real  part  of  a  complex  function,  U(t)  is  the  unit  step  function,  and  ap^ ,  Pk  y  Yp  ^ 
and  are  complex  constant  coefficients;  superscript  L  is  used  to  denote  the  linear  coefficient.  By  making  the 
proper  choices  of  complex  constant  coefficients  and  performing  Fourier  transforms,  we  can  readily  obtain  the 
familiar  Debye  and  Lorentz  forms  of  the  complex  conductivity  and  permittivity  in  the  frequency  domain.  [11] 

By  making  use  of  the  our  newly  formulated  three-dimensional  FDTD  algorithm  introduced  in  this  paper,  we 
investigate  much  needed  information  about  the  extent  of  penetration  of  microwaves  into  concrete  for  various  types 
of  externally  applied  incident  pulses.  We  discuss  the  results  of  numerical  simulations  of  high-power  microwaves 
penetrating  into  concrete  that  is  identified  by  frequency-dependent  conductivity  and  permittivity.  The  resulting 
FDTD  calculations  provide  second  order  accuracy  in  both  space  and  time  and  use  the  memory-efficient  recursive- 
convolution  (RC)  technique  to  update  electromagnetic-field  quantities.  We  present  in  this  paper  the  results  that  show 
the  distinct  difference  in  the  penetration  of  the  wide-band  pulse  as  compared  to  the  narrow-band  pulse  due  to  the 
dispersive  nature  of  concrete. 

n.  FDTD  FORMULATION 


In  light  of  Eqs.  (1.1)  and  (1.2),  we  write  Maxwell’s  equations  inside  the  dispersive  material  as 


d[fimt;x)] 


=  - Vx£(f;xJ, 


with 


K  K 

E(t;x)  +  E:^  ^P^p{t;x), 
p 

oo 

tp(t:x)^  \E(r;x)x<^>(t-T)dt, 


(2.1) 

(2.2) 

(2.3) 

(2.4) 


=  jE('i^;x)a^J’(t-T)dT,  (2.5) 

— oo 

where  is  the  magnetic  field  vector,  ju  is  the  magnetic  permeability,  Pp^(t;^  is  related  to  the  pth  terms  of  the 
linear  polarization  field  vector,  and  lx(t;x)  is  the  first-order  time-dependent  conductivity.  Using  an  FDTD 
algorithm,  the  above  equations  can  be  solved  numerically  at  each  time  step  provided  we  know  how  to  handle  Jk(U^ 
and  Pp(t;x)  numerically.  Therefore,  the  whole  solution  rests  on  the  question  of  how  to  carry  out  the  numerical 
evaluation  of  IkU;^  and  P_p(t;^  at  each  successive  time  step.  For  that  reason,  the  rest  of  this  section  is  devoted  to 
the  discussion  of  numerical  formulation  that  treats  and  Pp(t;x), 
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To  obtain  second-order  accuracy  in  time  in  evaluating  the  convolution  integrals,  E(t;x)  is  approxiniated  by  a 
piecewise  continuous  function  over  the  entire  temporal  integration  range  so  that  E(t;x)  changes  linearly  with  respect 
to  time  over  a  given  discrete  time  interval  [mAt,  (m+l)At],  where  m=0,l,...,n,  with  nAt  being  the  current  time  step 
[12-14].  Thus,  the  mathematical  expression  for  E(t;x)  takes  the  following  form  which  can  be  expressed  in  terms  of 
the  electric  field  values,  and  respectively,  evaluated  at  discrete  time  steps  t=mAt  and  t=(m+l)At  where 

these  two  successive  times  are  obtained  at  the  same  discrete  spatial  location  s={iAx,jAy,kAz)  with  Ax,  Ay  and  Az 
being  the  spatial  grid  sizes  in  the  jc,  y  and  z  directions,  respectively  (we  use  a  superscript  to  designate  the  discrete 
time  step  and  a  subscript  for  the  discrete  spatial  location): 


^  X  iAxjAy.kAi )  | 


F  + 

ft 


At 


for  0<mAt  <t  <(m-\rl )At  <( n+l )At 
for  t<0 


(2.6) 


For  three-dimensional  FDTD  calculations  in  Cartesian  coordinates,  we  need  to  solve  the  following  discrete  forms 
of  Maxwell’s  equations  that  obtained  from  Eqs.  (2.1)  and  (2.2)  by  finite  differencing  in  both  time  and  space  using 
the  usual  Yee  algorithm  [4]: 


)x-( ^y(k+'A>K~^ ^  ^ kJJ’ 

At 


At 

Ay 

At 


fi(  )y  -  fl(  ^  ^  ^  ^(i+^/2)3k  ^(i-^/2)jk  )z]  ^ij(k+y2)  )x  (  ^ij(k-‘/2)  k  ] I  ’ 


//(  -  //(  [  (  ^?Q^/2)k  )x]  ^(i+^/2)jk  )y  (  ^(i-^/2)jk  )y  IJ  ’ 

(DJ"  k  >.  =  k -<H’A  >y  -(<£«>,  ]) 


(2.7) 

(2.8) 

(2.9) 


-/  \  f’'^)x=(iAx,jAy.kAz)  K  j (  ^ k( £)x=^ iAx.jAy.kAz)  )x^  I  ’  (2  *0) 

Zii  K  K 

<DS;'k-(D;„),^t^[(H^\,).-(  >j~[<  HtXy:  >.  -( « "«»» 1) 

(n-H)At  (n+])At 

-{  J  ( E(t;x)^^i^,j^,tM)  )y  ]( t;x),^iax.jAy.ka,)  )ydt } ,  (2.11) 


(tt+l)At 


*■  nA 


( k  -<DI, ),  =  >,  -<  «("*«* 


(n+OAt  f"  1 

-{  \  (E(t;x)^^i^j^,tAz))zdt^(^'°^  j(JK(t;x),^iAx.jAy.kAz)>zdO-  (2-12) 

nAt  **  ^  nAt 

To  evaluate  the  right-hand  side  of  Eqs.  (2.10-2.12),  we  first  begin  by  substituting  Eq.  (1.3)  into  Eq.  (2.4)  then 
differentiating  these  integrals  with  respect  to  time  to  obtain  the  following  first-order  differential  equation  for 
complex  function  qJ' [14]: 


(n+])At 


3Gj(t;£) 

di 


+  l^Qp(  t;x  )=apE(f:x). 


(2.13) 


From  the  above  equation,  the  linear  polarization  vector,  Pj"  (t;x),  can  be  obtained  simply  by  taking  the  real  part  of 
Qp(t',x)-  Now,  solving  the  equation  exactly  in  time  by  using  the  integrating  factor,  exp(Ypt),  and  then  performing 
the  time  integration  from  nAt  to  (n+Z  )At,  we  have 
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njil+At 

n^+/it;x)=exp(- n/it;x)+apexp(- )  Jff  T;x)exp[-y^(  nAt-‘r)]dt .  (2.14) 

nAt 

Based  on  the  piecewise  linear  approximation  which  we  have  considered  for  the  temporal  behavior  of  we  can 
substitute  Eq.  (2.6)  into  the  right-hand  side  (i.e.,  the  inhomogeneous  part)  of  Eq.  (2.14)  for  E(TyX),  Evaluating  at 
spatial  location  ^^(iAxJAyykAz),  these  equations  result  in  the  following  recursive  relationship  for  (Qp^ in  terms 
of  (Qj^  )ijk\Ei/ and  Eyr': 

Q^^(nA+di;iA>c,jAy,kAzMQ‘'^)i:‘  ^exp(-Y^^At)(Q^^^,  ’Ip  (2-15) 

where 

( n+/ )At  At  L 

(¥pj))=ap  jexp[-y‘-p(nAt+At-x)]dz=ap^exp(-^pt)dT=^[l-exp(-r^pAt)] ,  (2.16) 

nAt  0 

-J-[i 

Y^At 


exp(-7iAt)]}  . 


(2.17) 


(¥p.i)=^l(At-T)exp(-/;,T)dT=^{] 

n  yp 


When  Eqs.  (2.15)  is  used  in  Eq.  (2.3)  to  obtain  discrete  forms  of  Eq.  (2.3)  at  successive  times  t-nAt  and 
/=(n+7jzdt  and  at  specific  spatial  location  ^{iAxyjAyykAz),  we  can  express  Dy/  and  Dy/"^^  ,  respectively,  in  the 
following  forms: 

ai=e.«. 

P 

and 

D^;’  =e,€^  Re{'^lexp(-r^dt  l^Sp(Vp.o  -EHjk  )(¥p., )1 } ■  (2-19) 

P 

Subtracting  Eq.  (2.18)  from  Eq.  (2.19),  we  obtain  the  following  expression  for  the  left-hand  side  of  Eqs.  (2.10-2.12) 
in  terms  of  (gp%"  ,  Eijk  and 


+  e,Re{Y  llexp(-Y^pdt )-!][( O'; ^  -M^k  X ¥p,  )1 }  •  (2-20) 

P 

To  evaluate  the  second  term  of  the  right-hand  side  of  Eqs.  (2.10-2.12),  we  substitute  Eq.  (2.6)  for 
E(t;2^xMiAx,)Ay,kAz)  Obtain  the  following  expression: 


(n+})At 

{  j  dtY^(T<J»  }  =  +E"j,  . 

nAt  ^ 


(2.21) 


To  evaluate  the  third  term  of  the  right-hand  side  of  Eqs.  (2.11-2.13),  we  first  substitute  Eq.  (1.4)  into  Eq.  (2.5) 
and  then  differentiate  the  integral  with  respect  to  time  to  obtain  the  following  first-order  differential  equation  for 
complex  function  /_/  (t;^: 


dlXiUx) 

dt 


+  e‘-l’;^(t;x)  =  fil;E(t;x). 


(2.22) 


As  in  the  case  of  the  linear  polarization  vector,  the  current  density  vector,  j/  (t;x)y  can  be  obtained  simply  by  taking 
the  real  part  of  L_k(U^  as  seen  in  the  above  equation.  Once  again,  we  solve  the  above  differential  equation  exactly 
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by  using  integrating  factor  exp(6^  t)  and  then  performing  the  time  integration  from  nAt  to  t  to  get  the  following 
expression  for  complex  function  l_J'  (t;x): 

t 

l';(t;x)  =  exp[e!^(nAt-t)]l‘;,(nAt;x)+fi!^exp(-el^t)^exp(e!;T)E(r;x)dr.  (2.23) 

nAt 

Again,  we  substitute  Eq.  (2.6)  into  the  above  equation  for  E(t;x)  based  on  the  piecewise  linear  approximation.  Then 
we  integrate  Eq.  (2.23)  from  nAt  to  nAt+At,  evaluated  at  spatial  location  ^(iAx,jAy,kAz),  to  obtain  the  following 
expression  for  the  third  term  of  the  right-hand  side  of  Eqs.  (2.10-2.12): 

(n+l)At  (n+l)At 

j Lk(  ^^  ^^x=(iAxJAy,kAz)^^  ^  j Lk( iAxJAy,kAz)  I 

nAt  nAt 

=  Re{ exp(eKnAt)l_\( nAt;iAx,  JAy.kAz)  J exp( -6^  t )dt 


( n+l  )At  t 


+  J^-Xp[0^(T  ^  )]  ]^ iAxJAy.kAz)  } 

nAt  nAt 

(n+J)At 

-Re{  expi  B^nAt  X  /  ^  )••*  |  exp(  -6^  1  )dt 


i  n+}  )At  t 


jdt  j  expl$^(T-l)][E^j^  + 


At 


=Re{J-[l-exp(-Bi;At)](t.  )^  (2-24) 

eb 


where 


in-i-l)At  t  M  1 

\dt]exp[d^(T-t)]dT=^{l-^[l-exp(-e^At)]},  (2-25) 

nAt  nAt  k  k 

^  Zi,  nAt  t'*'  ^ 

The  expression  for  (Lk  hj"  >  which  appears  in  Eq.  (2.24),  can  be  obtained  in  the  recursive  form  when  we  solve 
Eq.  (2.22)  exactly  for  J  J'(t:x)  using  integrating  factor  exp(0‘' t)  and  then  performing  the  time  integration  from 
t=(n-l)At  to  t=:nAt  while  making  use  of  the  piecewise  linear  approximation  for  E(t;x)  to  evaluate  the  “^tep-al  in  the 
right-hand  side  of  Eq.  (2.22).  The  result  is  that  we  obtain  the  following  recursive  relationship  for  (/^  )ijk  which  is 
expressed  in  terms  of  (7./  hjk  ' .  Eyk'^  and  : 

(i^)lk  ^expi-^tAXt.  <2.27) 


nAt  A  nL 

jexp[-Bt(nAt-r)]dr  =  fi^  jexp(-el^r  )dT’=t^[  l-exp(-e^At)J . 

(n-l)At  0  *■ 

( ^K.t  )^^\(dt-r}exp(  -e^n' )  dr= t-J—J  )]} . 

At  w  V*.  ^  JC 
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Finally,  when  Eqs.  (2.20),  (2.21)  and  (2.24)  are  substituted  into  Eqs.  (2.10-2.12),  we  obtain  the  following 
expressions  for  (Eijk*‘)x  ,  (Eijk"*^)y  and  (Eiji^*‘)z  in  terms  of  known  quantities  (Qp)ijk  ,  (L^kk"  >  Eijk"  and  Hijk'^'''' 
which  are  calculated  at  previous  time  steps  t=nAt  and  t=(n+}i)At: 

(2.30) 

(2.31) 

(2.32) 

where  aox,  0(^,0^,  and  ai  are  given  by 

+  Re{ Idt  )^k 

K  K 

-  (E^,  ),+e„  Ref  Ix'Z^expi- T^At)-!]  }+s„  (E^,  )xRe{Y,t(vU^-(Vp.>  ’  <^.33) 

P  P 

aoy  ^f^[(  Hftf^)jk  k  -( )x  -( )J}^^(E^k 
+  Ref  1(1^  }^'Y^^ii.exp{-ej^At)]}MEtjk  )y 

K  K 

-  €„e^(E^j,  )y+e.  Ref  [(Q%]yJ^[exp(-Y^pAt)-l]  (E^,  )yRef  ^[(VP^,o)-(¥p., )]  I  ^234) 

P  P 

aoz  (H^^A)k  )y  -(H^tl)jk  )y]}+^(El,  )f^CF<J» 

+  Ref  /(/^  ],Y,-^ll~^^P(-^K^t)]}+(Efjk  )z  EefY.K^K^M^';,,,)]} 

K  “k  K 

-  6„e„  (E^,  Ref  l(Q%]^^[exp(-Y^pAt)-l]  }+^„  (£^  KRef^[(Wpjo)-(Wp.,)]}>  (2-35) 

P  P 

a,  Ref'^(¥p.i)}-  (2-36) 

K  K  p 

In  our  FDTD  approach,  we  basically  need  to  update  Eqs.  (2.7-2.9),  (2.15),  (2.27)  and  (2.30--2.32),  respectively, 
for  Hyk"*'^  .  (Qp^)ijk"*‘ ,  (L^k"  and  at  each  time  step  in  order  to  carry  out  a  complete  three-dimensional 
computer  simulation  of  the  electric  field  response  in  linear  dispersive  materials  that  contain  both  fi-equency- 
dependent  electric  conductivity  and  polarization  terms. 

m.  SAMPLE  FDTD  CALCULATIONS:  ONE-DIMENSIONAL  CASE  STUDY  FOR  THE 
PENETRATION  OF  MICROWAVES  THROUGH  CONCRETE 

To  understand  dispersive  effects  in  penetration  through  concrete,  we  first  perform  one-dimensional  simulations 
of  narrow-band,  wide-band  and  ultrawide-band  electromagnetic  pulses  penetrating  through  finite-thickness  concrete 
using  the  newly  formulated  FDTD  algorithm.  We  perform  separate  simulations  for  narrow-band,  wide-band,  and 
ultrawide-band  incident  electromagnetic  pulses  to  investigate  the  effect  of  the  pulse  shape  in  the  transmitted  wave. 
For  these  simulations,  we  launch  a  time-dependent  electromagnetic  pulse  in  free-space  at  8.46  cm  (5  grid  cells)  from 
the  left-hand  edge  of  the  computational  space  and  propagate  the  pulse  from  left  to  right  in  the  positive  x-direction. 
After  traveling  16.83  meters  (995  grid  cells),  the  pulse  is  first  incident  on  20.30  cm  (12  grid  cells)  thick  dispersive 
concrete  that  is  characterized  by  first-order,  time-dependent  electric  susceptibility  and  conductivity.  After 
penetrating  through  dispersive  concrete,  the  wave  continues  to  propagate  in  Ifee-space  for  another  16.71  meters  (988 


(Elf^)y=-(aoy/aj), 

(E^;‘  ),=-(aod  a,). 
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grid  cells)  until  it  reaches  the  right  hand  edge  of  the  computational  space  that  is  terminated  by  a  LIAO  absorbing 
boundary  condition  [15].  Also,  the  LIAO  absorbing  boundary  condition  is  used  at  the  left-hand  edge  to  absOTb  all 
reflected  waves.  To  model  the  dispersive  behavior  of  concrete  at  microwave  frequencies,  we  consider  Debye 
forms  for  the  first-order,  time-dependent,  electric  susceptibility  and  conductivity  functions  [see  Eqs.  (1.3)  and  (1.4)] 
by  expressing  a/,  Pk  ,  and  in  terms  of  the  following  real  constants: 


=<5(€^ -e^)  and  7^  =  ^, 
=m )  and  6^  , 


(3.1) 

(3.2) 


where  <5 is  the  first-order  electric  susceptibility  damping  fi-equency,  e^is  relative  permittivity  at  ®  =  0,  e„is  relative 
permittivity  at  to  =  infinity,  t5  is  the  first-order  conductivity  damping  frequency,  <J,  is  electrical  conductivity  at 
(0  =  0,  and  o»is  electrical  conductivity  at  ©  =  infinity.  In  our  sample  FDTD  calculation,  we  treat  the  case  where  Ae 
dispersive  properties  of  concrete  are  represented  only  by  one  first-order,  time-dependent,  electric  susceptibility 
function  and  one  first-order,  time-dependent,  conductivity  function;  namely,  we  set  p„uix=  1  and  1  m  our 

FDTD  algorithm  obtained  in  the  previous  section.  Also,  we  set  <Tk to  be  equal  to  so  that  the  electnca 

conductivity  contribution  is  properly  taken  into  account  at  O)  =  0.  u  a\  u 

For  our  simulations,  we  investigate  three  different  incident  pulses:  the  sine-wave  pulse  (narrow-band),  the 
frequency-band  filtered  sine  function  pulse  (wide-band),  and  the  Gaussian  pulse  (ultrawide-band).  We  use  the 
following  mathematical  expressions  for  launching  of  the  three  incident  pulses: 


-  For  the  sine-wave  pulse:  Incident  Pulse  (t)= A  sin(a}J )[  U(  t  )-U( 


(3.3) 


-  For  the  frequency-band  filtered  sine  function  pulse: 

cc 

Incident  Pulse(t)=A  {  sine  [((0jt-t^^i^^)]-^sinc  }U(t),  (3.4) 

a 


-  For  the  Gaussian  pulse: 


Incident  Pulse  (t)-  A  exp[ 


(t-t 


delay 


^  width 


) 

-]U(t), 


(3.5) 


where  A  is  the  amplitude  of  the  electric  field  and  tjeiay  is  the  delay  time  for  the  incident  pulse.  We  arbitrarily  assign 
A  to  have  the  value  of  1.0  volt/meter,  to  have  values  of  4.23x10'’  second  (150  times  that  of  the  ume  increment) 
and  3.387x10'”  second  (12  times  that  of  the  time  increment),  respectively,  for  the  sine-wave  pulse  and  the  Gaussi^ 
pulse,  tjeim  to  have  the  value  of  2.82x10'®  second  for  both  the  frequency-band  filtered  sine  function  pulse  and  the 
Gaussian  pulse,  and  ©a  and  ©b  to  take  values  of  1.128xlO”  radians/second  and  4.512x10  radians/second, 
respectively.  Figurel,  Figure  2  and  Figure  3  show  the  time-domain  and  fi-equency-domain  plots  of  the  three  incident 
waves  we  use  in  our  simulations.  . 

We  select  the  total  number  of  simulation  cells  to  be  2000,  ranging  from  x  -  -1000  to  x  -  1000,  ''''t® 
ftee-space/concrete  interfaces  located  at  x  =  0  and  x=8.  We  launch  a  pulse  into  space  5  grid  cells  to  the  nght  of 
the  left-hand  edge  (at  x  =  -1000)  and  the  pulse  travels  in  the  positive  x-direction.  We  use  the  LIAO  absorbing 

boundary  condition  [15]  at  both  ends  of  the  computational  space. 

Following  are  the  basic  FDTD  parameters  that  we  use  in  our  one-dimensional  simulations; 


Uniform  grid  cell  size  (Ax)  =  1.691  centimeters. 

Total  simulation  distance  (Ax  times  the  total  number  of  cells)  =  33.82  meters. 

Time  step  increment  (At  =  Ax/2c)  =  28.2  picoseconds. 

Total  number  of  time  steps  =  5000  (total  simulation  time  =  0.141  microsecond ), 
where  c  is  the  speed  of  light.  . 

To  assign  values  for  parameters  that  describe  the  material  properties  of  dispersive  concrete  as  seen  in  bqs.  (J.i) 
and  (3.2),  we  make  use  of  the  concrete  data  shown  in  Figures  4  and  5  [1],  specifically  selecting  the  curves  at  the 
saturation  level  (S)  of  0.5. 

5.5;  and  S=  2.497x10'’ second;  =  0.05  mhos;  cr„  =  0.15  mhos;  15  =  2.0x10  second. 

Results  of  the  fully  dispersive  FDTD  calculations  are  shown  in  Figure  6,  Figure  7  and  Figure  8,  respectively,  for 
the  sine-wave  pulse,  the  frequency-band  filtered  sine  function  pulse,  and  the  Gaussian  pulse.  These  figures  are  snap 
shots  of  the  spatial  electric  field  patterns  taken  at  time  steps  of  1,500  and  3,500  (i.e.,  before  and  after  the  wave 
penetrates  through  concrete). 
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To  compare  results  of  our  fully  dispersive  FDTD  calculations  with  non-dispersive  FDTD  calculations,  we 
perform  nondispersive  FDTD  calculations  using  the  average  permittivity  and  conductivity  values  centered  on  the 
middle  frequency  of  the  incident  pulse.  Shown  in  Figure  9,  Figure  10  and  Figure  11  are  the  relative  cumulative 
energy  of  the  pulses  that  penetrate  through  concrete,  respectively,  for  the  sine-wave  pulse,  the  frequency-band 
filtered  sine  function  pulse,  and  the  Gaussian  pulse.  From  these  figures  we  can  see  that  the  wider  the  incident  pulse 
the  greater  the  discrepancy  between  fully  dispersive  and  non-dispersive  FDTD  calculations.  It  clearly  indicates  the 
importance  of  carrying  out  the  fully  dispersive  FDTD  calculations  for  dispersive  concrete  if  the  incident  pulse  is 
considered  to  be  relatively  broad  over  the  frequency  range  where  the  permittivity  and  conductivity  values  change 
significantly. 

IV.  CONCLUSIONS 

Based  on  the  FDTD  approach  we  presented  here,  we  can  solve  three-dimensional  Maxwell’s  equations  directly 
for  the  propagation  of  electromagnetic  waves  in  linear  dispersive  materials  that  exhibit  both  frequency-dependent 
electrical  conductivity  and  polarization.  Because  of  the  piecewise  linear  approximation  we  used  for  the  time 
dependent  part  of  the  electric  field  vector,  our  approach  provides  second  order  accuracy  in  time.  In  addition,  our 
approach  retains  all  the  advantages  of  the  usual  first-order  discrete  recursive  convolution  approach,  such  as  fast 
computational  speed  and  efficient  use  of  computer  memory; 

We  have  shown  that  it  is  critical  to  express  the  first-order,  time-dependent  electric  susceptibility  and  conductivity 
functions  in  the  exponential  forms  in  order  to  obtain  the  recursive  relationships  in  our  FDTD  algorithm.  An 
important  result  we  have  shown  in  this  paper  is  that  our  FDTD  algorithm  provides  second  order  accuracy  in  both 
time  and  space. 

Our  sample  one-dimensional  FDTD  calculations  show  that  it  is  very  important  to  include  fully  dispersive  effects 
in  the  analysis  of  microwaves  penetration  through  concrete,  especially  for  wide-band  and  ultrawide-band  incident 
pulses,  because  concrete  shows  significant  frequency-dependent  conductivity  and  permittivity  at  microwave 
frequencies.  Use  of  the  simple  frequency-averaged  conductivity  value  in  a  nondispersive  FDTD  calculation  will  not 
result  in  predicting  the  correct  amount  of  wave  energy  that  penetrates  through  concrete. 
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Electri  Field  (volt/meter) 


INCIDENT  WAVE  IN  THE  TIME  DOMAIN 


Figure  1:  Time-domain  (top  figure)  and  frequency-domain  (bottom  figure)  plots  of  a  sine- wave 
incident  pulse. 


10 


INCIDENT  WAVE  IN  THE  TIME  DOMAIN 
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Figure  3:  Time-domain  (top  figure)  and  frequency-domain  (bottom  figure)  plots  of  a  Gaussian 
incident  pulse. 
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Concrete  Real  Permittivity 


Figure  4:  Concrete  Permittivity  Data  [Reference  1] 

(note;  The  legend  in  the  left-hand  side  shows  the  specific  condition  that  the  data  is  evaluated 
for.  Also,  the  legend  in  the  right-hand  side  is  for  the  saturation  level  with  saturation 
levels  (S)  ranging  from  0. 1  (1 0%)  to  1 .0  (1 00%)  in  the  increment  of  0. 1 ) 
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Concrete  Conductivity  (mmhos) 
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Figure  5:  Concrete  Conductivity  Data  [Reference  1] 

(note:  The  legend  in  the  left-hand  side  shows  the  specific  condition  that  the  data  is  evaluated 
for.  Also,  the  legend  in  the  right-hand  side  is  for  the  saturation  level  with  saturation 
levels  (S)  ranging  from  0.1  (10%)  to  1.0  (100%)  in  the  increment  of  0.1) 
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Normalized  Electric  Field  Normalized  Electric  Field 


Thickness  of  concrete  is  1 2  grid  points  or  around  8  inches 


Figure  6:  Spatial  electric  field  profiles  for  a  sine-wave  pulse  taken  at  time  steps  =  1,500  (top  figure)  and 
3,500  (bottom  figure).  A  bold  line  in  the  middle  represents  the  thickness  of  concrete 
(12  grid  cells  or  8  inches). 


Normalized  Electric  Reid  Normalized  Electric  Reid 


Thickness  of  conaete  is  1 2  grid  points  or  around  8  inches 
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Figure  7:  Spatial  electric  field  profiles  for  a  fi-equency-band  filtered  sine  function  pulse  taken  at  time 
steps  =  1,500  (top  figure)  and  3,500  (bottom  figure).  A  bold  line  in  the  middle  represents  the 
thickness  of  concrete  (12  grid  cells  or  8  inches). 
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Normalized  Electric  Reid  Normalized  Electric  Reid 


Thickness  of  concrete  is  12  grid  points  or  around  8  inches 
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Figure  8:  Spatial  electric  field  profiles  for  a  Gaussian  pulse  taken  at  time  steps  =  1 ,500  (top  figure)  and 
3,500  (bottom  figure).  A  bold  line  in  the  middle  represents  the  thickness  of  concrete 
(12  grid  cells  or  8  inches). 


Transmitted  Energy/I nddent  Energy 


Transmitted  Energy/Incident  Energy 


OVERLAY  PLOTS  OF  NORMALIZED  TRANSMITTED  ENERGY 


Time  Steps  (Single  Time  Step=28.2  picoseconds) 


Figure  10:  Overlay  plots  of  cumulative  power  penetrated  through  concrete  for  fully  dispersive  FDTD 
and  nondispersive  (with  the  average  conductivity  value)  FDTD  calculations  in  the  case  of  a 
frequency-band  filtered  sine  function  pulse. 
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Transmitted  Energy/Incident  Energy 


OVERLAY  PLOTS  OF  NORMALIZED  TRANSMITTED  ENERGY 


Time  Steps  (Single  Time  Step=28.2  picoseconds) 


Figure  1 1 :  Overlay  plots  of  cumulative  power  penetrated  through  concrete  for  fully  dispersive  FDTD 

and  nondispersive  (with  the  average  conductivity  value)  FDTD  calculations  in  the  case  of  a 
Gaussian  pulse. 
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